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1 ABSTRACT 

O ■ 

High energy hadronic interactions can produce a final state char- 
Q_i| acterized by minijets separated by a large gap in the rapidity 

Qh! distribution of the produced secondary particles. We discuss the 

process by keeping into account the possibility of having mul- 
tiple parton collisions in the hadronic interaction. At Tevatron 
energy the correction to the single scattering term induced by 
the presence of multiparton interactions is large for transverse 
momenta smaller than 6GeV. 







1. Introduction 

In the kinematical regime of semihard hadronic collisions the momentum 
transfer q t is large enough to apply perturbation theory but it is kept fixed with 
the cm. energy, in such a way that the Regge limit t/s — > is reached from 
the perturbative side. The main feature, which becomes more and more apparent 
when approaching the Regge limit, is the increasing complexity of the process. In 
fact at large q t several major aspects of the interaction are described by a single 
partonic collision represented with a Feynman diagram at the lowest order in the 
coupling constant. When moving towards the semihard regime the Regge limit 
is approached both by the hadronic and by the typical partonic collision. As a 
consequence the partonic process is not well represented any more by means of a 
tree level Feynman diagram. Partonic interaction in the Regge limit have been 
investigated extensively [1], [2]. An approach which has received a lot of attention 
is the BFKL Pomeron[2]: the partonic reaction is described by the exchange of 
a gluon ladder, with vacuum quantum numbers in the t-channel, which is con- 
structed by neglecting the transverse momentum components when compared to 
the longitudinal ones. As a result of the simplified kinematics the steps of the 
ladder are ordered in rapidity and are build up with two basic elements: 
a- the gauge independent non-local vertices, which sum up the dominant term, 

in the t/s — > limit, of the diagrams with gluon emission from all near-by 

lines, and 

b- the Reggeization of the ^-channel gluons, which is introduced in order to keep 
into account the leading virtual corrections and which allows a solution to the 
infrared problem. 

The iteration of the ladder in the t-channel, is expressed as an integral equation 
(Lipatov's equation) which can be solved explicitly. By considering the parton 
process represented with a cut BFKL Pomeron one may write explicitly the ex- 
pression for the inclusive cross section to produce minijets, which are to be iden- 
tified with final state partons with transverse momentum larger than the lower 
threshold of observability q™ m . One may also consider the elastic parton inter- 
action with BFKL Pomeron exchange. The elastic parton collision gives rise to a 
distinct signature in the final state generated by the parton process: two minijets 
are produced with no further particles in the rapidity interval between them. On 
the contrary in the typical inelastic parton collision, represented with a cut BFKL 
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Pomeron, the gap is filled uniformly on the average by gluons. The dependence 
of the parton cross section on the rapidity gap is predicted by the perturbative 
calculation and it is related by unitarity to the dependence of the inelastic process 
on the width of the rapidity interval. The identification of these features in a 
semihard hadronic process would signal the underlying parton dynamics in a dis- 
tinctive way [3]. A difficulty which has been pointed out is that the hadronic event 
which contains the semihard partonic interaction is going to fill the gap in most 
of the cases. The perturbative cross section with rapidity gap has therefore been 
multiplied by a survival probability factor [4]. The survival probability has been 
recently estimated by taking explicitly into account the underlying soft hadronic 
event and it turns out to be roughly constant as a function of the rapidity gap [5]. 
The underlying event would therefore provide only a rescaling factor to the cross 
section calculated perturbatively and the behaviour of the actual experimental 
cross section as a function of the rapidity gap would still be linked directly to the 
BFKL dynamics. 

While testing the validity of the BFKL approach to semihard parton dynamics 
is presently one of the main topics in perturbative QCD, the delicate point is 
to keep properly into account the structure of the whole hadronic interaction, 
whose effect may mask the BFKL dynamics[6]. In fact the BFKL regime requires 
qmin sma ji w ith respect to all longitudinal momenta and smaller values of q™ in 
correspond to larger values for the partonic cross section. In the BFKL regime 
unitarity corrections are therefore important. Indeed the closer is parton dynamics 
to the BFKL limit the stronger is the effect of unitarization[6]. 

In the present paper we discuss the effect of the unitarization of the semihard 
hadronic interaction on the cross section for minijet production with rapidity gaps 
in the distribution of final state secondaries. The unitarization of the hadronic 
semihard interaction induces multiple semihard partonic collisions in the inelastic 
event and the cross section to produce minijets with associated rapidity gap is 
therefore modified. To keep into account multiple BFKL Pomeron exchanges we 
assume the validity of the AGK cutting rules for semihard interactions. As a 
consequence the whole semihard hadronic process is represented as a probabilistic 
superposition of multiple BFKL Pomeron exchanges. 

The paper is organized as follows: In the next paragraph the single parton in- 
teraction is discussed. The argument of the following section is multiple parton 
collisions; we describe the general formulation to the problem and we derive the 
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expression for the cross section with rapidity gap in the simplest case of multiple 
parton interactions. A numerical example and some concluding remarks are the 
argument of the last section while a few possible generalizations are discussed in 
the appendix. 



2. Single scattering term 



By considering the cut BFKL Pomeron one may write the expression for the 
cross section where two gluons interact producing many gluons and two of them, 
the ones with largest rapidity (in absolute value) in the overall cm. frame, are 
observed. If y is the separation in rapidity and k a , k^ are the transverse momenta 
of the observed gluons, the inclusive cross section is expressed as 



da j 



d 2 k a d 2 ki, 



CaOs 
U2 



f(k a ,k b ,y) 



(i) 



where Ca = N c is the number of colors, a s is the strong coupling constant and 
f(k a , ki,y) is the inverse Laplace transform of the solution to Lipatov's equation. 
Actually: 
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f(k a ,k b ,y)= V e in M due ^)y e iuin(kl/kl) 

[Z7Y) K a Kb n= _ oc J-00 

where 4> is the azimuth angle between the observed gluons, 



(2) 



u}(v,n) = —2 M 



n 



n I +1 



+ iv 



(3) 



and 



dlnT(z) 

dz 



(4) 



is the Digamma function. The inclusive cross section for production of two mini- 
jets, as a result of cutting the exchange in the forward direction of a BFKL 
Pomeron, is obtained by folding Eq.(l) with the structure functions of the in- 
teracting hadrons A and B: 



dx A dxtpk„<Pk h ~ feff(XA ' ffitefte*' k2b) d 2 td 2 h 



(5) 
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where f e f / is the effective structure function 



f eff (x) = G(x) + + QfW] ( 6 ) 

/ 

namely the gluon structure function plus 4/9 of the quark and anti-quark struc- 
ture functions with flavor /. In the BFKL dynamics one can relate the rapidity 
y of the minijet which carries most of the momentum of the initial state parton 
(primary minijet) with the fractional momentum variable x of the incoming par- 
ton. The relation is x = k t e y / 'y/s for forward final state partons with transverse 
momentum k t and rapidity y and x = k t e~ y / \fs for backward partons. Eq.(l) can 
be integrated on the transverse momenta down to the lower cut off q™ 171 : 



where 



with u){v) = 0) as a consequence of the integration on 0. &l is the inclusive 
cross section for minijet production in a parton process represented by a cut BFKL 
Pomeron. A different possibility which one may consider is to produce two minijets 
without cutting the BFKL Pomeron[7], namely by elastic scattering of two partons 
which exchange a BFKL Pomeron at momentum transfer larger than the lower cut 
off q™ 1 ™. The corresponding integrated partonic cross section is expressed as [3]: 

a S=(^) 4 1( ^F s (y) (9) 

and Fs(y) is the convolution of two BFKL propagators in the transverse momen- 
tum plane: 



/ min\2 

Fs(y) = / d 2 k 



16tt 3 



n miri 



J d 2 qd 2 q'f k (q,q',y) 



(10) 



k is the overall momentum exchanged through the ladder and the label S refers to 
the singlet exchange in the elastic partonic collision. Obviously one might consider 
also the possibility of having a octet exchange, the contribution to the cross section 
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is however sub leading at large rapidities [3]. The integrals on q and q' in Eq.(8) 
can be performed when using the integral representation of f k (q, q', y): 

J d*qd*qif k {q,q\y) = ±J dv ^ _^ /4)2 exp [ u (u)y] (11) 
and one obtains the asymptotic behavior at large y 

where z = a s CAy/^- For comparison the asymptotic behavior of &l is: 

(7 L (y) ~ — — — mm . 2 — — , /2 (13) 



tt y 2( ? r n ) 2 V(|c(3)^) 1/2 . 

The 'elastic' partonic cross section can be identified with the 'diffractive' 
cut of the double BFKL Pomeron exchange contribution to the forward parton 
amplitude[8]. At the same order, in the number of exchanged BFKL Pomerons, 
one needs therefore to keep into account also the one BFKL Pomeron cut and two 
BFKL Pomeron cut contributions to the partonic process. In fact the relative rate 
of elastic and inelastic parton processes, as expressed in Eq.(12) and in Eq.(13), is 
not consistent at large rapidity intervals because of the too rapid rise of the elastic 
cross section. We do not try to solve here this unitarity problem. To have an 
indication on the boundaries of the kinematical regime where the problem needs 
to be faced we have taken the simplest attitude. In analogy to the s-channel 
unitarization of the soft Pomeron exchange we have included in the semihard 
partonic interaction the exchange of two BFKL Pomerons and we have used the 
AGK cutting rules [9] to obtain the inelastic contributions to the cross section. The 
semihard partonic cross section duiy) is therefore expressed as 

vh{v) = °s{y) + (<Ti(y) - ±°s{y)) + 2a s (y) (14) 

where the single BFKL Pomeron exchange contributes with aL(y) and the con- 
tributions from the double BFKL Pomeron exchange, according with the AGK 
cutting rules, are: &s(y), the 'diffractive' contribution i.e. neither Pomeron is cut, 
— 4as(y), the one Pomeron cut, and +2as(y), the two Pomeron cut. Eq.(14) al- 
lows one to define the kinematical region of applicability of the approach. Indeed 
the one BFKL cut Pomeron contribution to the cross section must be positive[10]: 
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(a L (y) - 4a s (y)) > 



(15) 



Eq.(15) fixes a limiting value to the rapidity interval y as a function of as, as 
an example for a s = .18 and as = -22 one obtains about 14 and 10 rapidity 
units respectively. Strictly speaking the contribution expressed by Eq.(15) implies 
that the one-cut-Pomeron events and the two-cut- Pomeron events are distinguish- 
able in term of their multiplicity. If the fluctuations in multiplicity prevents this 
identification we get the milder condition 



which requires only the distinction between diffractive and non-diffractive events. 

We therefore express the semihard cross section <Jh{v), corresponding to a 
single partonic interaction, as 



and CFp{y) is the contribution from production of secondaries, both from one and 
two cut BFKL Pomerons. 

3. General Framework for Multiparton Interactions 

Given the scale q™ in one may separate parton fluctuations with a lifetime long 
with respect to q™ in and parton fluctuations with a lifetime short with respect to 
g ( mm . Parton fluctuations which have a long lifetime are of non-perturbative origin 
and can be associated to the initial state. Parton fluctuations with a short lifetime 
can be treated perturbatively and can be associated to the semihard interaction. In 
the simplest case, when the typical value of the rapidity interval y in the partonic 
interaction is not too large, the semihard partonic interaction is described at the 
lowest order in the number of exchanged BFKL Pomerons. If q™ m is relatively 
small one faces however a unitarity problem even if one is in a regime where the 
partonic interaction is well described within perturbation theory because of the 
large value of the integrated semihard hadronic cross sectionfll]: by integrating 
Eq.(5) with the cut off q™ in one obtains a cross section which easily exceeds the 
value of the total cross section. On the other hand the partonic cross section, 



(& L (y)-2& s (y)) >0 



(15') 



oh{v) = 6s{y) + {6 L {y) - 2a s (y)) = a s (y) + a P (y) 



(16) 
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Eq.(16), is in comparison still rather small. The large value of the integrated 
inclusive cross section is therefore the consequence of the large flux of partons 
in the initial state, which gives rise to an average number of partonic collisions 
larger than one[12]. In the typical semihard hadronic process different partonic 
interactions are localized in the transverse plane in different regions, of size of order 
{&h) inside the overlap volume of the two interacting hadrons, whose transverse 
size is of a few fm 2 . At Tevatron energy, considering partonic interactions with 
cross section a H , as expressed in Eq.(16), and q™ in = 6GeV : the average distance 
in rapidity between produced minijets is about five units. The corresponding value 
of the partonic cross section is ch{{v)) — 8 x 10 _2 m6 and one may correspondingly 
estimate that in a central collision the average number of partonic interactions is ~ 
2^4. These features characterize the kinematical regime which we are presently 
interested in: actually the 'elementary' partonic interaction is well described within 
the BFKL dynamics, in such a way that the bound in Eq.(15) is satisfied, and the 
inclusive hadronic minijet cross section is larger than the inelastic cross section, 
in such a way that the rate of multiple parton interactions is sizeable. 

To discuss multiple parton collisions we follow the approach used in ref. [6] and 
[13,14]: In the case of soft interactions multi-Reggeon exchanges are conveniently 
taken into account by making use of the AGK cutting rules [9]. Although no gen- 
eral proof of their validity is available in the case of semi-hard interactions, it has 
nevertheless been possible to show that the cutting rules hold for one of the com- 
ponents of the interaction which is leading in the large-s fixed-t limit[15]. If one 
assumes the validity of the cutting rules for semi-hard interactions, one is allowed 
to represent the semi-hard cross section ou as a probabilistic distribution of mul- 
tiple semi-hard parton collisions [14]. The most general expression for an requires 
however the introduction of the whole infinite set of multiparton distributions[16], 
which keep into account hadron fluctuations in the parton number. To that pur- 
pose we introduce the exclusive /c-body parton distribution W"( fc ) (u\ . . . Uk) , namely 
the probabilities to find a hadron in a fluctuation with k partons with coordinates 
ui...Uf~, Ui = (bi,Xi) standing for the transverse partonic coordinate (pi) and 
longitudinal fractional momentum (xi). From the exclusive distributions and the 
auxiliary functions J{u) one may construct the generating functional Z[J}: 

Z[J]=^ — I J(u 1 )...J(u n )W n (u 1 ,...u n )du 1 ...du n (17) 
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in such a way that the exclusive distributions are the coefficients of the expansion 
of the generating functional around J = 0, while the inclusive distributions are 
the coefficients of the expansion around J = 1. A very general expression for the 
semihard cross section, which is consistent with the cutting rules, is therefore: 



Here the (3 is the impact parameter between the two interacting hadrons and 
&H(ui,u'j), represents the probability for the parton i of the A-hadron to have a 
semihard interaction with the parton j of the S-hadron. The semi-hard cross sec- 
tion is constructed by summing over all possible partonic configurations of the two 
interacting hadrons (the sums over n and m) and, for each configuration with n 
A-partons and m S-partons, summing over all possible multiple partonic interac- 
tions. This last sum is constructed by asking for the probability of no interaction 
between the two configurations (actually nlLiIl^iil ~~ a i,j] )• The difference 
from one of the probability of no interaction gives the sum over all semi-hard 
interactions. o~h((3) is then the probability to have at least one semihard parton 
interaction when the impact parameter in the hadronic collision is equal to (3. The 
semi- hard cross section is obtained by integrating the probability o~h(I3) on the 
impact parameter. Analogously, the elementary semi-hard cross section &h(x,x') 
is obtained by integrating the elementary interaction probability 0"#(w, u') on the 
relative transverse coordinate b — b'. 

In Eq.(18) oh is constructed by summing all possible semihard two-body 
parton collisions. Multiple semi- hard parton collisions are of two distinct kinds, 
disconnected collisions and rescatterings. In a disconnected collisions different 
pairs of partons interact independently at different points in the transverse plane. 
In a rescattering a high energy parton interacts several times, with momentum 
exchange larger than q™ 171 , with different target partons and all interactions are 





(18) 



n 



m 




i=lj=l 
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localized in the same region, approximately of size (l/q , J nm ) 2 , in the transverse 
plane. The most important contribution to the semi-hard cross section arises 
from the disconnected partonic collisions. In fact, at a given number of partonic 
collisions, the incoming parton flux is maximized in the configuration where all 
collisions are disconnected. It might therefore be meaningful to obtain a simpler 
expression for o H by neglecting the rescattering processes in Eq.(18). To that 
purpose we expand the interaction probability ( the factor in curly brackets ) as 
sums and suppress all addenda containing repeated indices: 

{l-n[l-*y]}=}>5](7y--5] Wkl + ... (19) 

i,j ij ' ij k&,tyj 

Because of the symmetry of the derivative operators in Eq.(18) one can replace 
the expression in Eq.(19) with: 

nm&n — 2\ n ( n ~ ^)m(rn — l)<3"ii<5"22 + • • • (20) 

in such a way that the sums over n and m in Eq.(18) can be performed explicitly. 
As a consequence the cross section at fixed impact parameter &h{0) is expressed 
by the operatorial form: 

* H (p)= [i-exp(-— ■ a ■ -)]z A [J + 1}Z B [I + 1} \ j=i=o (21) 

where the dependence on the variables u and u' is understood. 

The expression of cth(P), as given by Eq.(21), is still too complicated to 
be worked out, since all possible multi-parton correlations are implicitly present 
in Z. The simplest possibility is to neglect all correlations in the multi-parton 
distributions. In this case one writes 



Z[J+1] = exp J D(u)J(u)du 



(22) 



where D(u) is the average number of partons. The cross section assumes therefore 
the eikonal form: 



a H = J rf 2 /?[l-exp(-$(/3))] (23) 



where 



9 



rVM rVM 

m3) = $s(P) + $ P (P)= / dy dy'(<l>s(P;y,y') + MP;y,y')) (24) 

J V m J V 



with 

4sAfcy,y') = J d 2 bD A {b,x{ y ))a s Ay' '-y)D B (b- p,x'( y ')) (25) 

and umi Vm are the maximum and minimum rapidity values allowed by kinemat- 
ics. The index S in Eq.(24) refers to the elastic parton interaction, the singlet 
exchange in Eq.(16), and the index P to the inelastic interaction, the one and 
two cut Pomeron contributions in Eq.(16). The different contributions from mul- 
tiple semihard parton collisions to the cross section are explicit if one expands the 
exponential in Eq.(23) as follows: 

x dyxdy'JfeiPiyiiy'x) + (f>p(P;yi,y[)) , na . 

x / dy v dy v [(\) S {^y v ,y l v )+(\)p{^y v ,y v ) 
Jy v <v' v V 

One is interested in the component of an which represents two minijets at 
rapidities y and y', in the central rapidity region, with associated gap Ay = y' — y 
in the rapidity distribution of secondary produced gluons. To that purpose one 
needs to exclude in Eq.(26) both the 'elastic' terms, with final state minijets in 
the gap, and all the inelastic partonic interactions, generated with elementary 
probability (pp. In fact the cut BFKL Pomeron originates soft gluons distributed 
on the average uniformly in rapidity. Local fluctuations in the rapidity distribution 
of soft gluons, which could leave the gap empty also in the case of an inelastic 
elementary parton interaction, may be neglected if one considers a reasonably 
large rapidity gap. The cross section to observe two minijets at rapidities y and 
y', with the gap Ay = y' — y in the rapidity distribution of secondaries, is therefore 
expressed as 
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da H (Ay) 
dydy' 



where 



d 2 (3 



2_^v<ps{f3;y,y)- j — - 



00 ry ryM 

+ / dy<t> s {P;y,y') / dy'<t> s (p;y,y' 



x 



[MA Ay)] 



i/-2 



77! 



* s (J3,Ay)= / dy / dy' 'Mfrvd) 
•Jym Jy' 

After summing on v one obtains 



-*(£) 



(27) 



(28) 



da H (A y) 
dydy 



y 



MPW,y')+ I d y( ps(P;y,y') dy'M^y^'n 

Vm Jy' J (29) 



x 



exp{$ s (/3; Ay) - $s(0) - 



The two addenda in Eq.(29) are the single and double 'elastic' scattering 
terms. In the single scattering term both observed minijets are produced in the 
same elementary partonic interaction, in the double scattering term the two mini- 
jets are generated in different partonic collisions. Both terms are multiplied by the 
absorption factor exp{ — ($>s(P) — ®s(P'i Ay))} that removes the 'elastic' parton 
interactions which would fill the gap, actually those which produce minijets with 
rapidities y and y' such that y < y or y' < y'. At a fixed value of f3 the cross 
section is multiplied by exp{— $p(/3)} which is the probability of not having any 
inelastic partonic interaction in the process. One may recognize in Eq.(29) the 
semihard contribution to the survival probability factor (S 2 (f3)} of ref.[4]. Actu- 
ally exp{— &s(/3) — ®p(/3)) } is the probability factor of not having any semihard 
activity in the underlying hadronic event. A more detailed analysis of the origin 
of the survival factor and of the suppression factor due to some elastic scattering 
is presented in the Appendix, in the context of a multiparton dynamics. In the 
same Appendix the influence of the possible deviation from the Poissonian form 
of the original parton distribution is also discussed. 
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4. Discussion and conclusions 

High energy hadronic interactions, with production of minijets and associated 
rapidity gap in the distribution of secondaries, is a process where one would expect 
to observe explicitly the unitarity relation at the level of partonic amplitudes. A 
closer inspection shows that the unitarity problem which one faces is twofold. 
On one hand, if one believes that the semihard interaction between partons can 
be described with the BFKL approach, one needs to solve the problem of the 
too rapid rise of the partonic cross section as a function of the rapidity interval 
between primary minijets. On the other hand, even if the size of the partonic 
cross section is not too large, namely if the hadronic cm. energy is such that the 
typical rapidity interval between primary minijets is relatively small, the semihard 
hadronic cross section may still be too large. One can therefore distinguish three 
different regimes: 

I- The cutoff is sizeable with respect to the typical energy available to the semi- 
hard partonic interaction. The corresponding 'elementary' parton interaction 
is small, no unitarization is needed and the semihard cross section is well 
described by a single partonic collision. 
II- The cutoff is moved towards relatively smaller values, with respect to the 
available energy. A single partonic interaction is still well described by the 
BFKL dynamics. The semihard hadronic cross section is however too large 
with respect to the total inelastic cross section and unitarity corrections are 
to be taken into account. The unitarization of the hadronic semihard cross 
section is achieved by taking into account multiparton interactions, namely 
different pairs of partons interacting independently with BFKL Pomeron ex- 
change. Typically the different partonic interactions are localized at different 
points in the transverse plane, in the region of overlap of the matter distribu- 
tion of the two hadrons. 
Ill- With even smaller values of the cutoff one may still be in the regime where 
perturbative QCD can be used, since the value of ots{qT in ) is small, but 
the typical rapidity interval between primary minijets is too large so that 
the 'elementary' parton process is not well described any more by the single 
BFKL Pomeron exchange and also the 'elementary' parton process has to be 
unitarized. One may obtain an indication on the limits between regions II 
and III by testing whether the bound in Eq.(15) is satisfied. 
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To have a quantitative indication on the boundaries of the kinematical regions, we 
have worked out a numerical example. Since the present available information on 
multiple parton interactions is limited to the scale factor which gives the rate of 
double parton interactions [17], while no information is available on the multipar- 
ton correlations, we have considered the simplest possibility, namely the Poisson 
distribution for the multiparton distributions. On the grounds that the main con- 
tribution to the multiple parton interactions is represented by the disconnected 
partonic collisions, we have neglected parton rescatterings. The unitarized expres- 
sion for the cross section is therefore given explicitly as a function of the input 
which is used to evaluate the single scattering term, namely the average number 
of partons D(b, x) and the 'elementary' partonic cross section an- We have factor- 
ized D{b,x) as / e //(x) x F(b), where f e ff(x) is the effective structure function as 
expressed in Eq.(6) and F(b) is a Gaussian, normalized to one and such as to give 
for the double scattering term the scale factor a e ff (see Eq.(26) of Ref.[6]) con- 
sistent with the experimental indication[17]. In our numerical example we have 
chosen a e ff = 20mb and as a scale factor for the structure functions we have 
taken q™ m /2. as is a free parameter in the BFKL approach, one expects however 
that the value of as which one should use is not too different from the value of 
the running as at the scale of the typical momentum transferred in the process. 
We have chosen as a value of as the value of the running coupling computed with 
qmin 1 2 as a sca l e factor. The values of the semihard cross section an, as expressed 
in Eq.(23), which we obtain with this input are consistent with the experimen- 
tal values published by UA1[18]. At each value of the hadronic cm. energy the 
boundary of the kinematical regions I, II and III are identified by the choice of 
the cutoff q™ m . The curve which corresponds to larger values of m fig-1 nas 
been drawn requiring that the unitarized hadronic semihard cross section an, as 
expressed in Eq.(23), is 20% smaller with respect to the single scattering term. 
The lower curve corresponds to the value of q™™ which, on the average, namely 
after integration with the structure functions, saturate the bound in Eq.(15). The 
two curves identify the three regions mentioned above. 

Moving from large values of q t to the semihard region one faces therefore 
two different unitarity problems, which signal the appearance of different levels 
of structure in the hadronic interaction. The gap in the rapidity distribution of 
produced secondaries is an effect which derives from the unitarity relation applied 
to the 'elementary' parton amplitude. The regime where 'elementary' interactions 
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with rapidity gap are a sizeable component of the 'elementary' parton process is 
however a regime where the overall hadronic process is already structured in a 
non-trivial way, because of the large amount of multiparton interactions. To have 
a quantitative feeling of the effect of multiparton interactions, we have performed 
a numerical calculation in the simplest example already considered to obtain the 
curves in fig.l. The expression of the cross section as a function of the rapidity 
interval Ay is given in Eq.(29). We have considered pp interactions at Tevatron 
energy and as a lower threshold to observe minijets we have taken the value q™ 171 = 
5GeV. In fig. 2 we plot the cross section, as expressed in Eq.(29), divided by 
the survival probability factor exp{— $>s(P) — $p (/#)}• The continuous curve is 
obtained by using as a input the value cr e // = 20m6 and the dashed curve is 
obtained by using a e ff = Ylrnh. The dotted curve is the contribution of the single 
scattering term alone. 

As it is shown in fig. 2 the effect of unitarization on the behaviour of the cross 
section is large. In the actual case the main modification to the dependence on Ay 
is due to the presence of multiple 'elastic' parton scatterings whose effect on the 
cross section is twofold. A different dependence on Ay, with respect to the single 
scattering term, is induced by the presence of two different sources. The first is the 
contribution of the process where the two observed minijets originate in different 
'elastic' partonic interactions, the second term in Eq.(29), the second source for 
the different dependence on Ay is the correction induced by multiple 'elastic' 
scatterings to the survival probability factor. In fact not all underlying hadron 
activity needs to be excluded. 'Elastic' parton scatterings which produce minijets 
outside the gap are allowed and the corresponding contribution to the cross section 
depends on Ay. The effect of the inelastic semihard partonic interactions is, on 
the contrary, factorized at fixed impact parameter (3 and independent on Ay. 
The main effect of the inelastic partonic processes is to contribute to the survival 
probability (S 2 ) of ref.[4] rather than modifying the dependence on Ay. 

The region where q™ in is relatively small is also the region which is closer to 
the BFKL kinematics and where, as a consequence, the BFKL approach to parton 
dynamics is better justified. In high energy hadronic interactions, at relatively 
small values of q™ m , unitarity corrections are however large and have to be taken 
into account. The physical effect which underlies the need of unitarization is the 
increasing complexity of the semihard interaction which appears when q™ %n is low- 
ered, as it is shown in fig.l. In the kinematical region which has been discussed in 
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the present paper, corresponding with region II in fig.l, the 'elementary' parton 
process is well described within the BFKL approach. The typical hadronic interac- 
tion is however characterized by several 'elementary' partonic collisions, localized 
at different points in the transverse plane. One of the reasons of interest in the 
actual kinematical regime is that, in region II, the non perturbative component of 
the process, which is factorized in the multiparton distributions, represents a new 
piece of information on the hadron structure with respect to the hadron structure 
functions of large pt physics. 

Appendix 

A.l General features and inelastic scattering 

In this Appendix the way in which an overall term independent of Ay and 
further corrections explicitly dependent on Ay may arise out of the many parton 
dynamics is discussed more in detail. The functional formulation, which has been 
already employed, is used again and more systematically. The longitudinal mo- 
menta are always given in terms of the fractional momentum x, the relation with 
the corresponding rapidity y was mentioned in Sec. 2. 

If we look to configurations where there are two radiated partons (jets) with 
rapidity gap sitting in the central rapidity region, in the two-hadron cm., then 
the two partons belong to different original hadrons. The parton population is 
divided into three segments x a ,Xb,x c : the x a interval correspond to the partons 
which would fill the gap, so to the parton that, in the configuration looked for, are 
not scattered; the Xb correspond to the partons outside the gap, so to the partons 
which are simply not observed; the x c correspond to the parton observed at one 
end of the gap; the segment in which x c lies will be ideally shrunk to one point. 
The three regions of variations of x reflect into three regions of variations of u, 
but clearly the impact parameter b is not affected by this operation. 

This decomposition suggests a way of rewriting the generating functional 
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which is sometimes more convenient : 

Z[J] = y ] — [ J{u\) ■ ■ ■J(u n )W n (u 1 ,...u n )du 1 ...du n 

111 f 

= Y] -t-ttt / J(u a )---J(u b )---J(u c )Wp !q!k (u a ,...u b ,...u c )du a ...du b ...du c , 

(AA) 

p is the number of partons of kind a, q is the number of partons of kind 6, k is the 
number of partons of kind c, and an analogous expression is written for Z[I(u')]. 

We start considering the effect of inelastic processes and we express the prob- 
ability that: 

I- A parton u c scatters elastically against a parton u' c or a parton u c scatters 
elastically against a parton u' h and parton u' c scatters elastically against a parton 
Ub thus giving rise to the two partons at the end of the gap and nothing is produced 
in the middle. 

II- No inelastic scattering takes place, since such a process would give rise 
also to partons which could fill the gap. 

In formulae the requirement I is expressed by an operator T while the re- 
quirement II is expressed by an operator both applied to the product Z[J]Z[I]. 
The actual form of the two operators is: 



du c du' c 



8 c {u)a E {u c , u' c )S' c (u')+ 



5 c (u)a E (u C: u' b )8' h {u') 5 b (u)a E (ub, u' c )5' c (u')du b du h 



(A2) 



n,n' 




-^6 1 (u)---8 n (u)---8' 1 (u')---8' n ,(u')x 
n\ n'\ 



(A3) 



jQjl — bi n {u, u')}dudu 



In the definition of Q there is no point in distinguishing the intervals of ra- 
pidity. Some notational simplification is obtained by setting: 



5J(u r 



= 8 m (u) 



81{u' m - (3) 
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= 8' m (u') m = a, 6. 



(AA) 



A first step in the choice of some specific distribution could be to neglect the 
genuine many-body correlations, in this case the generating functional is in fact 
reduced to an ordinary function of a linear functional of the sources 



z[j] = <f>(y) y = J D(u)j(u)du 



(A.5) 



The division of the field of variation of x into three parts induces a corresponding 
decomposition y = y a + y b + y c . 

Even with the restriction to a pure one-body density the general expression is 
complicated. A very relevant and effective simplification is produced if we assume 
a Poissonian distribution for the initial partonic distributions 



z[j] =ex P [y-y ] 



(A6). 



The normalization term y o = j D(u)du ensures that Z[l] = 1. 

In looking for some general features of the inelastic rescattering it is however 
possible to use a general partonic distribution without many-body correlations, 
as it is described by the generating functional $>(y) ; so we let the functional 
differential operators Q and T act on the product $(3^ J]) $(3^-7]) at the end the 
auxiliary sources are put to zero: J = I = 0. From the action of the operator Q 
one obtains the intermediate result: 



n,n' 




1 1 



n\ n'\ 



A)- D(u n ) ■■■D(u' 1 - /?)••• D(u' n , - (3) 



x 



- a In (u,u')}dudu'& n \y[J])& n '\y[I}) 



(A.7) 



From the further action of the operator JF and the condition J = I = one 
obtains the final result: 



K (P)=J2 I I ^D(u 1 )...D(u n )...D(u' 1 -(3)...D(u' n ,-(3) 

n,n' 

JJjl _ &i n (u,u')]dudu' X 

'D(u c )a E (u Cl u' c )D(u' c -/3)^ n+1) (0)^ n ' +1) (0)+ (A8) 

D(u c )a E (u c , u' b )D(u' h - 13) D(u b )a E (u b , u' c )D(u' c - (3)du b du b 
$ (n+2) (0)$ (n ' +2) (0)l 



/ 



du c du' c 
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So in general the expression has the form 



K(P) = • + G^W) ■ F»m (A.9) 

n,n' 

For a Poissonian partonic distribution it results for every derivative $( n )(0) = 
exp[— y o ], so the double sum over n, n' can in principle be carried out yielding an 
overall factor G{j3) = J2 n ,n' [^n,n' (P) + ^n[n'0^)] w hich multiplies the remaining 
term 



F(p)= J du c du' c [D{u c )a E {u c ,u c )D(u' c - (3)+ 

J D(u c )a E (u c , u' b )D(u b - 13) D(u b )a E (u b , u c )D(u c - f3)du b du' b 

(A10) 

In other words it results that, at fixed hadronic impact parameter (3 the inelastic 
processes give simply rise to a multiplicative factor to the fundamental amplitude. 

In presence of another kind of distribution the treatment is less straightfor- 
ward because the sum over n, n' does not allows the extraction of a common term 
F((3). It seems however likely that in the expression of K the second addendum, 
containing F H is more important than the first one, because it involves a less 
exclusive condition, if this is true then the sum 

n,n' 

can be still be carried out yielding again a multiplicative factor in front of the 
fundamental amplitude at fixed hadronic impact parameter. 

This essential simplicity is destroyed if we are in presence of sizeable two- 
body correlations. To be definite we may consider an example where we have a 
Poissonian distribution corrected by a two-body correlation 

z[j] = exp[y-y ]. 

y = J D(u)J{u)du + \j C(u,u)J{u)J{u)dudu. (All) 

Then in performing the derivatives in order to calculate K{j3) we end unavoidably 
with expressions where the correlation term C(u, u) links the inelastic component 
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with the elastic fundamental one and so the decomposition given in Eq.(A.lO) 
is no longer possible, in this case we expect that the inelastic processes modify 
strongly not only the size but also the shape of the fundamental elastic process. 



A. 2 Effect of the multiple elastic scattering 

It is clear that in presence of a large flux of partons the inelastic processes 
that have been considered till now are not the only disturbing effects i.e. there are 
other dynamical processes which also act in the sense of filling the rapidity gap 
that the "fundamental" process described in the present formalism by T produces. 
With reference to the notation introduced in eq (A.l), one must take into account 
also the following elastic processes: 

A parton u a which scatters against a parton u' a , this event will fill the gap. A 
parton u a which scatters against a parton u' b or a parton u b which scatters against 
a parton u' a , this event will also fill the gap. The statement that these kind of 
collisions must not happen is expressed by a veto operator which, in absence of 
inelastic interactions would have the following form: 

£ = / / ^^M...S b (u)...S' a (u')...S' b (u')x 

p,p' q,q' J J y ' H ' F ' H ' 

ED 1 ~~ ^E{u a , u' a )] ■ [1 - & E {u a , %)] • [1 - & E (u b , u' a )]du a du b du' a du' h 

(A 12) 

The fact that the parton of kind b are unobserved suggests that the sum over 
g, q' is the first operation to be performed. 

For a fixed p the term containing the sum over q' is 

A(tt a> i, • • • , u a , P ) =5 a (ui) ■ ■ ■ 5 a (u p ) x 

E^n[/ t 1 - & E(u a ,l,u' k )]---[l -& E (Ua,p,u' k )}S' b (u' k )du' k . ^ AA ^ 
q> q ' k J 

Since all the u' k are dummy variables the sum over q' can be carried out and gives 
an exponential; then the product of the [1 — a] factors in the exponent may be 
expanded in terms containing no a, one a, two a, and so on. 
A(u 0) i, • • • , u a , p ) = 8 a {ui) ■ ■ ■ 5 a (u p ) x 

J { l ~Y^ a E{u a j, u' k ) + \a E {u a , u u' k )a E (u aJ , u' k ) + • • ^8' h {u' k )du' k , 

(A14) 



exp 
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The sum over q is obviously treated in the same way and we get for 8 the expres- 
sion: 

-j -77 Yli 1 - <^K, u' a )]A(u aA , • • • , u a!P )A(u' a l , • • • , u' a pl ) Y[du a du' a 

p, P ' 

(A15) 

In order to obtain more explicit expressions it is necessary, now, to expand both 
the operators A and the term — a] in multiple interactions, but at this step 
it seems that the treatment becomes clearer if one chooses some definite form of 
the parton distributions. 

We choose now a Poissonian distribution of partons * and let the functional 
differential operator S act on the product Z[J]Z[I], after the action of T . Only at 
the end the auxiliary sources are put to zero: J = I = 0. As it was already said, 
for the operator £ one is forced to proceed in steps corresponding to successive 
reinteractions of the same partons: from the operator A one gets: 



A(u a , 



D(u a ,i) ■ ••D(u a ,p)x 



exp 



f^-Yl °E(u a ,i, u' b ) + \a E {u a ,n u' h )a E {u a ^, u' b ) H ]D{u' h - (3)du' b 

i i^j 

(A16) 

Also the term P = ^[1 — &E(u a ,i, u ' a , r )] =Yl[^~ &i,r] is expanded with the result: 

P = 1 - &l > r + Yl Yl h & i,rVj,s "EE EtVSO^r^^M + ' ' ' • (^-17) 

In the repeated sums the pairs of indices cannot be equal e.g. (i, r) 7^ (j, s); it is 
however possible to have either i = j or r = s, what correspond to a rescattering 
of a particular parton. If we took only the term 1 and the simple sum in the 
exponent of A, out of the operator S we would get this result: 



E =exp 
exp 



D(u a )F(u a - [5)du a 



exp 



D(u' a - (3)T(u' a )du' a 



x 



D(u b )du b + / D(u' b )du' b 



(A18) 



* In order to perform the actual calculation it may be useful to remember the 
identity f(d/dx)g(x)\ x= o = g(d/dx)f(x)\ x =o, which is evident whenever the two 
functions admit a power expansion around 0; in particular it yields f(d/dx)e hx \ x =o - 
f(h) 
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In this formula the absorption factors V have been introduced, in their definition 
the property that a depends only on the difference of the impact parameters has 
been used. 

T(u a - j3) = exp - J a E (u a , u' b )D(u' b - /3)du' b ^ 



T(u' a ) = esxp - J a E {u' a ,u b )D{u b )du' b ^. 



The factor E Q takes into account the interaction of the partons of kind a with 
those of kind 6, but not the interaction of the a-partons among themselves. It 
could be put in a form that could be of easier interpretation, i.e. 



E Q = exp 



exp 



x 



D(u)du + J D{u')du' 
- [ D(u a )[l-T(u a -f3)]du a 



exp[-y D(u' a -p)[l-T(u' a )]du' a 

(A18 ; ) 

In this expression the first exponential, where also the small contribution from 
u c has been included, takes simply away the normalization term exp[— y o — y f ], 
it does not belong to the multiple scattering, while the second term, E' Q is an 
"absorption" term due to some of the disconnected collision which fill the gap. We 
can do better and consider the other terms in the expansion of P = — o"i,j], 
with the restriction of excluding the repeated interaction of the same a-parton, in 
formulae i 7^ j and r 7^ s in this way a further factor is produced which multiply 
the former expression, i.e. 

■ [ D(u a )T(u a - p)& E {ua,u' a )D{u' a - (3)T(u' a )du a du' a } (A19) 



Ei = exp 

Going on means to produce expressions quite elaborated, one could, as an example, 
calculate the factors which describe a double interaction of partons either of kind 
a or of kind b. When the partons a are allowed to collide twice, but not three 
times, the following factor is produced: 

E 2 = exp \ jD(u a )T(u a - /3)a E (u a , u' a )D{u' a - 0) x 

T(u' a )a E (u' a , Ua)D(ua)T(ua - (3)du a duadu' a + 

(A20) 

D (u a - P)T{u a )aE\u a , u a )D{u a ) x 
T(u a - P)a E (u a , u' a )D(u' a - (3)T(u' d )du a du' a du'^ . 
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The expansion can be further continued, e.g. by considering also the collision 
of two partons of kind b and then the various results can be collected. From 
the term explicitly calculated one gets: E ps E' EiE 2 ■ ■ ■. This factor modifies 
qualitatively the result expressed in Eq.(A.9): because in that expression the terms 
G n , n ' were independent of the rapidities defining the gap which appear in F n>n / 
through u c = (x c , b); now, on the contrary these variables appear explicitly in the 
terms E. 

Other form of partonic distributions give rise to more complicated final expres- 
sions, but some of them allow anyhow some elaboration, in particular as already 
remarked in [14], the generating functional of a negative binomial distribution 
allow a representation which is an integral transform of a Poisson generating func- 
tional: 

So, in principle, from the previous expressions holding for a Poisson distribution 
of partons it would be possible to get the analogous one for a negative binomial 
distribution by applying the following prescriptions: 

I- Take away the normalization term, A/p = exp[— y o — y' ] and substitute it 

by 

^ B = [r(a)]-Hi-y r(i-y r 

II- Multiply by a factor t every D(u) explicit or implicit, i.e. inside Y and B; 
and multiply by a factor s every D(u') explicit or implicit, i.e. inside Y and 

B. 

III- Multiply the overall resulting expression by e~ t ~ s (ts) a ~ 1 . 

IV- Integrate the result in dtds from to +00. 

It is evident that the Poissonian and the negative binomial distribution for 
the original partons are only the simplest ones which can be chosen; in term of 
integral transforms of the generating functions other distributions could be built 
up. 

A. 3 Elastic and inelastic scattering 

The real case implies evidently both elastic and inelastic scattering. The 
requirement that there be no production process filling the gap is expressed by an 
operator of the form 
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Yl[l-aT(u a ,u' a )]-[l-a T (u a ,u' b )}-[l-a T (u b ,u' a )} 
Y\[t - vin(ub, u' b )]du a dubdu' a du' h 

a T = a E + cfin (A. 23) 

The statement expressed by this operator is that while the parton of kind a must 
not suffer interactions at all, the partons of kind b may interact provided their 
interaction is elastic. In its full form this expression appears quite intractable, 
however one can follow the idea that the most relevant dynamical feature is the 
inelastic scattering and the elastic is a perturbation. It does not seems convenient, 
however, to start from the inelastic operator Q, it is more efficient to start with 
an "exaggerated" form of the operator C such that also the elastic interactions of 
the partons b are excluded, by defining 

n,n' 

(A24) 

which looks very similar to Q, eq(A.4) and so yields, by itself, a factor independent 
of y. Then one has to correct the result by an expansion in a E that uses eq (A. 23). 
Since the coordinates u are all equivalent the expansion has the same combinatorial 
structure as: 

[A + B] qq ' = A qq ' + qq'BA^-W-V . A q ~ x ■ A q ~ x + B 2 ■ ■ ■ 

In this way, through a shift in q and q', the integrand of the operator C a is repro- 
duced with additional factors. The actual form of the first order correction in a E 
is: 

p,p' q,q' 

Y[[l-a T (ua,u' a )] ■ [1 - a T (u a ,u b )] ■ [1 - a T (u b ,u' a )] ■ [1 - a T (u b , u' b )\ x 
du a du b du a du b 

& E {vb, v' b ) JJ[1 - &T(ub, v' b )] - & T (vb, u' b )]5b(v)5' b (v')dv b dv' b . 

(A25) 
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In a limit in which terms like &e ■ &t are neglected, while keeping the powers of 
&t the correction may be factorized in the form 

Ci « C ■ J a E {v b , v' h )8 b {v)5' h {v')dv b dv' b . (A.26') 

This approximation for the correction C\ suggests also a generalization. The ex- 
pression itself originates from the division of the whole flux of parton in a pair 
suffering only elastic scattering and in a rest for which all interaction are consid- 
ered. This attitude can be extended to more general situation, on the ground that, 
as already noticed, in presence of large parton fluxes the multiple disconnected in- 
teractions may become more relevant than the reinteractions. So, with reference 
to Eq.(A.22) we can divide the q partons of kind b into two subsets of r and s 
partons, obviously with r + s = q, and the same for q' . The coordinates of the r 
partons will be still denoted by u, the coordinates of the s partons will be called 
v. In this way Eq.(A.22) takes the following form: 



/l 1 1 1 1 1 
v\7\ 7\V\ V\ V\ 5ai - U) ' " 5>b ^ ' lit 1 " <r]du a du b dv b du' a du' b dv' b , 



p,p' r,r' s,s 



(A.27) 

the product — a] contains nine kinds of factors. Now we decide, on the 
basis of the physical consideration outlined before that the v parton suffer only 
elastic scattering, without any other interaction, then necessarily s = s', and there 
are si ways of pairing the v with the v'. So it is not difficult to see that the rest 
of C reconstruct the a factor C a of Eq.(15), whereas the sum over s yields an 
exponential. In this way the approximate expression for C is given by: 

C w C ■ exp [J 8 b (v')a E (v b , v' h )8' h {v')dv b dv' b ^ . (A.2T) 

When the parton distribution is Poissonian out of Eq.(A.27') one gets the 
expression of the first correction due to multiple scattering depending on Ay which 
corresponds to what was previously stated in section 3, Eq.(29). Anyhow, in a way 
wholly independent of the validity of the approximation we see that the correction 
is always positive because it represents a compensation for the previous exclusion 
of configurations that should be allowed, the amount of the correction diminishes 
by enlarging the rapidity gap since the integration runs over the complementary 
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domain. For our purposes the most relevant feature of the elastic interactions is 
that they produce effects depending from the very beginning on the y of the gap 
because the integrations over v, v' depend explicitly on x c , x' c which are precisely 
the variables that set the difference between partons of kind a and partons of kind 
b. 
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Figure captions 



Fig.l- The three different kinematical regions which characterize semihard hadronic 
interactions. I: only the single partonic collision, described by a single BFKL 
Pomeron exchange, is relevant; II: multiparton collisions are to be taken into 
account, each partonic interaction is however well described by a single BFKL 
Pomeron exchange; III: the single BFKL Pomeron exchange is not any more 
an adequate description of the single parton interaction. 

Fig. 2- Cross section for production of minijets with rapidity g function of the 

gap Ay. The process is pp at yfs = 1.8TeV and q^ %n = 5GeV. The dotted 
curve is the single scattering term. The continuous and the dashed curves 
include all disconnected multiple parton collisions. Input to the continuous 
curve is a e ff = 20mb and to the dashed curve is cr e ff = 12m6 (see text). 
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